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ABSTRACT 


First-principles  density  functional  theory  (DFT)  calculations  using  the  generalized  gradient  approximation 
(GGA)  have  been  performed  to  study  the  adsorption  of  a  series  of  all-nitrogen  and  high-nitrogen  compounds  of 
increasing  sizes  and  complexity  on  the  Al(l  11)  surface.  The  calculations  employ  periodic  slab  models  with  4  A1 
layers,  ranging  in  size  from  (3x3)  to  (6x6)  surface  unit  cells,  and  containing  up  to  144  A1  atoms.  Complementary 
quantum  chemical  calculations,  utilizing  DFT  and  second-order  perturbation  theory  methods,  of  the  ground  state 
potential  energy  surfaces  of  the  corresponding  polynitrogen/high  nitrogen  species  in  the  absence  of  the  aluminum 
surface  also  have  been  performed.  The  initial  set  of  studies  performed  in  the  first  year  of  this  challenge  project, 
which  focused  on  the  adsorption  and  reaction  properties  of  Nx  (x=l-5),  NHX  (x=  1  -3),  N2HX  (x=l-4),  N3H,  N3H3  and 
N4H4  species  on  Al(l  11),  have  been  extended  to  include  larger  polynitrogen  systems  such  as  Nf„  N8,  Ni0  and  Ni2. 
The  C48N12  fullerene-like  molecule  was  included  in  this  series  in  place  of  N60,  which  was  found  to  be  unstable  in  the 
gas  phase.  Furthermore,  simple  nitrogen-containing  heterocycles  such  as  furazan  (C2N2H2O),  2H-1,2,3-  and  4H- 
1,2, 4-triazole  (C2N3H3),  and  1H-  and  2H-tetrazole  (CN4H2)  have  been  studied.  Additional  compounds  of  interest 
include  the  monosubstituted  di(cyclopropyl)triazole  (C| iN(,Hi2)  and  di(cyclobutyl)triazole  (Ci3N5Hi6)  derivatives  of 
triazene.  For  these  systems,  surface  interaction  mechanisms  involving  both  dissociative  and  nondissociative 
processes  have  been  characterized  as  a  function  of  molecular  orientation  and  surface  site.  The  dissociative 
mechanisms  generally  include  elimination  of  one  or  more  N2  molecules. 

I.  Introduction 

Investigation  of  the  interactions  of  energetic  high-nitrogen  compounds  with  ultrafine  or  nanophase  aluminum 
particles  is  a  topic  of  current  DoD  interest  with  potential  applications  in  rocket  and  missile  propulsion.  One  of  the 
topics  of  prime  importance  is  to  determine  if  the  thermal  instability  inherent  in  many  high  nitrogen  compounds  might 
be  reduced  by  adsorption  or  chemisorption  onto  the  surface  of  aluminum,  while  simultaneously  forming  a  protective 
coating  on  the  metal  surface  which  minimizes  the  formation  of  an  inert  oxide  surface  layer. 

II.  Computational  Method 

The  calculations  performed  in  this  study  were  done  using  the  Vienna  ah  initio  simulation  package 
(V  ASP). 13  This  program  evaluates  the  total  energy  of  periodically  repeating  geometries  based  on  density -functional 
theory  and  the  pseudopotential  approximation.  In  this  case  the  electron-ion  interaction  is  described  by  fully  non¬ 
local  optimized  ultrasoft  pseudopotentials  similar  to  those  introduced  by  Vanderbilt.4  5  Periodic  boundary  conditions 
are  used,  with  the  one-electron  pseudo-orbitals  expanded  over  a  plane-wave  basis  set.  All  calculations  have  been 
performed  using  a  cutoff  energy  of  435  eV  except  for  the  furazan  calculations.  In  this  case,  the  energy  cutoff  was 
increased  to  495  eV,  due  to  the  presence  of  O  atom  which  requires  a  higher  cutoff  energy.  These  cutoff  energies  are 
consistent  to  the  precision  level  high  required  by  the  pseudopotentials  used  by  VASP  code.  Calculations  were 
performed  using  the  generalized  gradient  approximation  (GGA)  density  functional  theory  with  PW91  exchange- 
correlation  functional.6  The  sampling  of  the  Brillouin  zone  was  performed  using  a  Monkhorst-Pack  scheme.7  The 
minimum  energy  path  between  different  minima  was  determined  using  the  climbing  image  nudged  elastic  band 
method  developed  by  Jonsson  and  co-workers.8'9 

The  GAMESS10  quantum  chemistry  code  was  used  to  predict  the  structures  and  properties  of  isolated 
polynitrogen  and  high-nitrogen  compounds.  The  theoretical  methods  employed  included  second  order  perturbation 
theory11  (MP2,  also  known  as  MBPT(2))  and  density  functional  theory  (DFT)  using  the  hybrid  B3LYP12  functional. 
The  6-31  lG(2df,p)  basis  set13  was  used  throughout. 
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III.  Results  and  Discussion 

In  our  systematic  analysis  of  high-nitrogen  and  polynitrogen  compounds  with  the  aluminum  surface,  we 
have  considered  a  diverse  number  of  chemical  compounds  in  order  to  determine  the  adsorption  and  reaction 
properties  of  different  chemical  functionalities.  The  first  set  of  compounds  considered  is  the  sequence  of  all-nitrogen 
molecules  N6,  N8,  N10  and  N12,  which  extends  the  previous  investigations  of  Nx  (x=l-5)  from  the  first  year  of  this 
project.  The  C48N12  fullerene-like  molecule  was  included  in  this  series  in  place  of  N60,  which  was  found  to  be 
unstable  in  the  gas  phase.  The  second  set  of  compounds  considered  in  this  study  is  comprised  of  simple  nitrogen- 
containing  heterocycles,  including  furazan  (C2N2H2O),  2H- 1,2, 3-triazole  and  4H-  1,2,4-triazole  (C2N3H3),  and  1H- 
and  2H-tetrazole  (CN4H2).  The  third  set  includes  the  monosubstituted  di(cyclopropyl)triazole  and 
di(cyclobutyl)triazole  derivatives  of  triazene  (CnN6H12  and  C13N5H16,  respectively.) 

Due  to  the  diversity  in  the  size  of  the  molecular  systems  analyzed  in  this  study  we  have  used  in 
computations  slab  models  of  different  sizes.  These  range  from  (3x3)  surface  supercells  used  in  the  case  of  small 
systems  to  (6x6)  surface  units  for  the  largest  systems.  In  all  cases  the  slab  models  contained  four  layers  of  atoms. 
The  optimizations  have  been  performed  for  both  the  adsorbate  and  the  A1  atoms  in  the  top  two  layers  while  the 
bottom  two  layers  were  frozen  at  the  bulk  optimized  conditions.  Adsorption  at  different  surface  sites  and  for 
different  molecular  orientations  have  been  investigated  and  the  corresponding  adsorption  energies  have  been 
calculated  using  the  expression  Eads  =  Emoiec  +  ES|ab  -  E(moiec+S|ab)  where  Emoiec  is  the  energy  of  the  isolated  molecular 
species  at  its  equilibrium  geometry,  Es|ilb  is  the  total  energy  of  the  isolated  aluminum  slab,  and  E(moiec+siab)  is  the  total 
energy  of  the  adsorbate/slab  system.  A  positive  Eads  corresponds  to  a  thermodynamically  stable  adsorbate/slab 
system.  In  the  present  study,  only  neutral  species  have  been  considered. 

Nx  (x=6,  8,  10,  12)  Adsorption  on  Al(lll).  The  adsorption  properties  of  Nx  (x=6,  8,  10,  12)  have  been 
characterized  using  a  (3x3)  supercell  aluminum  slab  model  for  N6  and  N8,  and  a  (4x4)  slab  model  for  N10  and  N12. 
Representative  adsorption  configurations  are  presented  in  Figures  la-11.  In  the  case  of  N6,  both  an  open-chain  (C2h) 
and  a  cyclic  (D2)  isomer  have  been  considered.  Similarly,  for  N8  both  the  azidopentazole  (Cs)  and  octaazapentalene 
(D2h)  isomers  have  been  investigated.  Upon  optimization,  the  open-chain  N6  isomer  was  found  to  dissociate  to  form 
two  N3  species  which  bind  on  the  surface  in  either  horizontal  or  vertical  orientations,  as  illustrated  in  Figure  la.  In 
contrast,  the  cyclic  N6  isomer  molecularly  adsorbs  on  the  surface  with  formation  of  multiple  Al-N  bonds,  as  shown  in 
Figure  lb.  For  cyclic  N6,  the  binding  energies  range  from  49.3  to  64.8  kcal/mol,  relative  to  the  isolated  slab  and  gas 
phase  N6  molecule.  In  the  case  of  N8,  both  the  azidopentazole  (Figures  lc-le)  and  octaazapentalene  (Figure  If)  were 
found  to  adsorb  molecularly.  For  azidopentazole,  significant  variations  in  binding  energies,  ranging  from  10.5  to 
63.1  kcal/mol,  were  found  as  a  function  of  specific  molecular  orientation  and  the  number  of  Al-N  bonds  formed  (see 
Figures  lc-le.)  In  some  instances,  azidopentazole  was  found  to  dissociate  on  the  surface,  leading  to  formation  of 
adsorbed  N5+N3  or  N6+N2  fragment  pairs.  All  of  these  reactive  processes  are  exothermic  with  respect  to  the 
undissociated  configurations  (see  Figures  lc- Id.)  Fikewise,  dissociation  of  octaazapentalene  to  N5  +  N3  (Figure  If) 
is  exothermic  by  about  64  kcal/mol.  Similar  reaction  processes  were  also  observed  for  the  N10  molecule.  For 
example,  N10  can  adsorb  without  undergoing  fragmentation  (Figure  lg)  or  by  a  dissociative  process  in  which  one  of 
the  molecular  rings  is  broken  with  subsequent  formation  of  multiple  Al-N  bonds  (Figure  lh.)  Additional  dissociative 
mechanisms  include  the  formation  of  two  adsorbed  N5  species  (Figure  li)  or  an  adsorbed  N6  fragment  plus  two  N2 
molecules  which  readily  desorb  (Figure  lj).  The  dissociative  adsorption  mechanisms  of  N10  were  found  to  be  highly 
exothermic,  with  dissociation  enthalpies  as  high  as  101  kcal/mol  relative  to  the  adsorbed,  undissociated  state. 
Finally,  in  the  case  of  the  N12  system,  molecular  adsorption  was  found  to  be  possible  with  binding  energies  in  the 
range  55.5-60.4  kcal/mol  (see  Figures  lk-11.)  The  presence  of  several  Al-N  bonds  in  this  case  contributes  to  the 
increased  stability  of  this  molecular  system. 

By  comparing  the  current  set  of  results  to  the  prior  adsorption  data  for  Nx  (x=l,5),  it  can  be  concluded  that 
adsorbed  N  atoms  are  the  most  stable  species.  Additionally,  for  the  larger  systems  analyzed  here,  the  surface 
dissociation  leading  to  formation  of  N,  N2,  N3  fragments  is  highly  exothermic.  Not  surprisingly,  the  interaction  of  N2 
with  the  A1  surface  is  weak,  leading  to  facile  desorption. 

Adsorption  of  C48Ni2  Fullerene-like  Systems.  Three  distinct  configurations  of  the  C48Ni2  fullerene-like  cage  on 
the  A1  surface  have  been  identified.  The  most  stable  structure,  with  a  binding  energy  of  12  kcal/mol,  is  obtained  by 
the  attachment  of  the  C48NJ2  cage  to  the  A1  surface  via  adjacent  carbon  atoms,  as  shown  in  Figure  In.  The 
configuration  in  which  the  fullerene  cage  is  attached  to  the  A1  surface  via  a  pair  of  neighboring  carbon  and  nitrogen 
atoms  (Figure  lo)  is  somewhat  less  stable,  with  a  binding  energy  of  10.9  kcal/mol.  The  least  stable  structure,  shown 
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in  Figure  lp,  involves  binding  of  the  fullerene  cage  to  the  surface  via  two  non-adjacent  carbon  atoms,  with  a  binding 
energy  of  4.5  kcal/mol.  No  dissociative  adsorption  pathways  were  observed  for  C48N12  on  the  A1  surface. 

Furazan  Adsorption  and  Dissociation  on  Al(lll).  Adsorption  of  the  furazan  ring  on  the  A1  surface  can  occur  in 
vertical  or  horizontal  configurations  as  illustrated  in  Figures  2a-2d,  with  a  preference  for  the  vertical  state  (Figure 
2a),  which  has  a  modest  binding  energy  of  7.3  kcal/mol.  The  adsorbed  furazan  ring  can  readily  dissociate  on  the 
surface  with  a  large  amount  of  energy  released.  This  is  illustrated  in  Figures  2m-2n  which  show  the  minimum  energy 
pathways  for  molecular  dissociation  on  the  surface  starting  from  an  initial  horizontal  local  minimum  (Figure  2m)  and 
from  an  initial  vertical  configuration  (Figure  2n.)  As  indicated  in  these  figures,  the  corresponding  activation  energies 
for  molecular  dissociation  are  very  small  but  the  overall  reactions  are  highly  exothermic;  e.g.,  187  kcal/mol  starting 
from  the  horizontal  configuration. 

Adsorption  of  2H-1, 2, 3-triazole,  4H-1, 2, 4-triazole,  lH-tetrazole,  and  2H-tetrazole.  Illustrative  examples  of  the 
adsorption  configurations  of  2H-1, 2, 3-triazole,  4H-  1,2,4-triazole,  lH-tetrazole,  and  2H-tetrazole  are  indicated  in 
Figures  2e-21.  The  corresponding  binding  energies  are  moderate,  with  values  of  10.8,  15.7,  10.7,  and  8.0  kcal/mol, 
respectively.  In  all  cases,  a  vertical  orientation  of  the  ring  and  attachment  to  the  surface  via  one  or  two  nitrogen 
atoms  is  the  preferred  adsorption  configuration. 

Adsorption  of  mono-substituted  triazenes.  As  an  extension  of  our  previous  results  on  the  aluminum  surface 
adsorption  of  the  1,3,5-  and  1,2,3-triazene  molecules,  the  surface  interactions  of  the  monosubstituted 
di(cyclobutyl)triazole  (Ci3Hi6N6)  and  di(cyclopropyl)triazole  (ChHi2N6)  derivatives  of  triazene  have  been 
investigated.  These  calculations  have  been  performed  using  a  (5x5)  surface  slab  model.  Figure  3  illustrates  the 
results  for  these  two  systems,  which  shows  four  adsorption  configurations  for  the  cyclobutyl  (Figures  3a-3d)  and 
cyclopropyl  (Figures  3e-3g)  derivatives.  In  all  cases,  the  cycloalkyl  groups  do  not  bind  to  the  aluminum  surface. 
The  most  stable  configurations  of  the  cyclobutyl  (Figure  3a)  and  cyclopropyl  (Figure  3e)  systems  involve  binding  to 
the  surface  via  the  formation  of  multiple  bonds  between  aluminum  and  the  6-membered  triazene  moiety  plus  an 
additional  Al-N  bond  to  a  nitrogen  atom  of  the  triazole  substituent,  with  binding  energies  of  25.1  and  21.6  kcal.mol, 
respectively.  A  similar  but  less  stable  configuration  for  the  cyclobutyl  derivative  was  found  in  which  there  is  no 
bond  between  an  aluminum  atom  and  the  triazole  moiety  (Figure  3b.)  Additional  configurations  in  which  bonding  to 
the  surface  occurs  via  a  single  nitrogen  atom  of  the  triazene  ring  are  also  possible,  as  shown  in  Figures  3d  and  3f, 
with  respective  binding  energies  of  12.3  and  12.5  kcal/mol.  The  least  stable  configurations  show  a  horizontal 
orientation  of  the  triazene  moiety,  as  illustrated  in  Figures  3c,  3g,  and  3h.  The  latter  two  structures  show  that  a 
significant  portion  of  the  stabilization  energy  gained  from  formation  of  Al-N  or  Al-C  bonds  can  be  lost  in  the 
deformation  of  the  gas  phase  structure  and/or  in  the  upward  movement  of  some  of  the  surface  A1  atoms.  This  is  the 
case  for  the  structures  depicted  in  Figures  3g  and  3h,  which  both  have  negative  binding  energies  despite  the  fact  that 
a  number  of  Al-N  and  Al-C  bonds  have  been  formed.  In  the  present  sign  convention,  the  negative  binding  energies 
indicate  that  these  structures  are  thermodynamically  unstable  relative  to  the  isolated  gas  phase  molecule  and  the 
isolated  surface. 

IV.  Conclusions 

A  systematic  study  of  a  series  of  polynitrogen  and  high  nitrogen  species  (Nx;  x=6,8,10,12,  the  C48N12 
fullerene  analogue,  simple  nitrogen-containing  heterocycles  such  as  furazan,  triazole,  and  tetrazole,  and 
monosubstituted  di(cycloalkyl)triazole  derivatives  of  triazene)  interacting  with  the  aluminum  (111)  surface  have 
been  investigated  using  plane-wave  DFT  calculations  in  conjunction  with  aluminum  surface  slab  models.  Our 
calculations  indicate  that,  in  the  majority  of  cases  molecular  adsorption  takes  place  by  formation  of  multiple  Al-N 
or  Al-C  bonds.  For  the  ensemble  of  molecular  systems  considered  here,  a  significant  range  of  variation  in 
binding  energies  has  been  determined  depending  on  the  number  of  surface  bonds  involved  as  well  as  the  surface 
site  and  the  molecular  orientation.  Additionally,  the  adsorption  was  found  to  lead  to  important  molecular 
deformations  and  to  upward  shifts  of  some  of  the  surface  aluminum  atoms.  For  the  all-nitrogen  species, 
dissociative  adsorption  readily  occurs,  particularly  for  the  N6-Ni0  species,  with  a  large  energy  release. 
Dissociation  often  leads  to  desorption  of  N2  species  and/or  deposition  of  N  atoms.  The  overall  dissociation 
processes  were  found  to  be  exothermic.  Such  findings  were  also  found  for  furazan,  where  molecular  dissociation 
takes  place  with  a  large  amount  of  energy  release. 

Future  work  will  extend  this  set  of  investigations  to  include  additional  high-nitrogen  and  polynitrogen 
systems  such  as  3,3’-azobis(6-amino-l,2,4,5-tetrazine)  (C4N12H4)  as  well  as  a  systematic  analysis  of  the  thermal 
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effects  on  the  adsorption  and  reaction  processes. 
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Figure  Captions 

Figure  1.  Side  and  top  views  of  the  most  stable  adsorption  configurations  of  (a)  Nfl  (open-chain)  and  (b)  N6  (cyclic) 
molecules.  Side  views  of  N8-azidopentazole  (c)-(e),  N8-octaazapentalene  (f),  Ni0  (g)-(j),  and  Ni2  (k)-(l)  on  the 
Al(l  11)  surface.  Panel  (m)  illustrates  the  gas  phase  structure  of  the  C48Ni2  molecule.  Adsorption  configurations 
involving  adjacent  carbon  atoms,  adjacent  carbon  and  nitrogen  atoms,  or  nonadjacent  carbon  atoms  in  C48Ni2  are 
shown  in  panels  (n)-(p),  respectively. 

Figure  2.  Side  views  of  representative  binding  configurations  of  furazan  (a-d),  2H-1, 2, 3-triazole  (e,f),  4H-1, 2,4- 
triazole  (g,h),  lH-tetrazole  (i,j),  and  2H-tetrazole  (k,l).  Panels  (m)  and  (n)  indicate  the  minimum  energy  pathways 
for  furazan  dissociation  starting  from  horizontal  and  vertical  adsorption  configurations,  respectively. 

Figure  3.  Representative  adsorption  configurations  of  the  monosubstituted  di(cyclobutyl)triazole  (a-d)  and 
di(cyclopropyl)triazole  (e-h)  triazene  derivatives  and  the  corresponding  binding  energies. 
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